ANL-HEP-PR-06-54 
EFI-06-09 

Gravitons and Dark Matter in Universal Extra 

Dimensions 



o 

^ ! Nausheen R. Shah" "^ and Carlos E.M. Wagner^ 
D 

I Enrico Fermi Institute"' and Kavli Institute for Cosmological Physics^, 

O '• University of Chicago, Chicago, IL 60637, USA 



OO 
O 

o 



^HEP Division, Argonne National Laboratory, 9700 Cass Ave., Argonne, IL 60439, USA 

February 2, 2008 



Abstract 



Models of Universal Extra Dimensions (UED) at the TeV scale lead to 
^ i the presence of Kaluza Klein (KK) excitations of the ordinary fermions and 

5h ' bosons of the Standard Model that may be observed at hadron and lepton 

^ colliders. A conserved discrete symmetry, KK-parity, ensures the stability of 

the lightest KK particle (LKP), which, if neutral, becomes a good dark matter 
particle. It has been recently shown that for a certain range of masses of the 
LKP a relic density consistent with the experimentally observed one may be 
obtained. These works, however, ignore the impact of KK graviton production 
at early times. Whether the is the LKP or not, the G" tower thus produced 
can decay to the LKP, and depending on the reheating temperature, may lead 
to a modification of the relic density. In this article, we show that this effect 
may lead to a relevant modification of the range of KK masses consistent 
with the observed relic density. Additionally, if evidence for UED is observed 
experimentally, we find a stringent upper limit on the reheating temperature 
depending on the mass of the LKP observed. 



1 Introduction 



Theoretical models constructed to understand physics beyond the standard model 
(SM) (most notably string theory) frequently imply the existence of extra dimen- 
sions. It has been a particular challenge in phenomenology to understand how these 
extra dimensions would be realized and manifested in our observable 3 + 1 dimen- 
sional world. The number, shape and size of these dimensions, as well as the particles 
allowed to propagate in them give rise to several different models, all having different 
phenomenological implications. 

We will be considering universal extra dimensions (UED) where the SM fields 
propagate in all the extra dimensions. For toroidal compactification, this would 
imply a tower of Kaluza-Klein (KK) particles for every SM particle, each carrying 
KK number. Momentum conservation in the extra dimensions implies KK number 
conservation. However, the requirement of obtaining the proper SM chiral modes 
at low energies leads to constraints on the possible compactification geometries. An 
example for D = 5 is the orbifold /Z2 where the Z2 projects out half of the zero 
modes, leaving only SM fields. Additionally, this breaks translation invariance along 
the extra dimensions, so KK number is no longer conserved. A residual symmetry, 
KK-parity, is still present and it is sufficient to ensure the stability of the lightest 
KK particle (LKP). It also ensures that KK particles are always produced in pairs, 
allowing for good agreement between theory and experiment for small values of the 
compactification scale, of the order of a few hundred GeV 

The stability of the LKP allows for an interesting candidate for dark matter 
(DM) The LKP is expected to be weakly interacting and electrically neutral if 
it is to be considered a candidate for dark matter. The usual candidate is , the KK 
partner of the hypercharge gauge boson. For mxK ~ O (1) TeV, -B^ gives excellent 
agreement with the observed relic density. There have been many further analyses of 
the relic density in UED, including a more proper treatment of coannihilation effects, 
and the impact of the inclusion of second KK level resonances jHI El 13 El El HPj • 

However, most of these studies ignore the gravitational sector. We shall mostly 
work within the context of one extra dimension. The gravitons couple extremely 
weakly and therefore are generally considered to be unimportant for collider studies. 
Cosmo logically, however, they can have a significant effect ^, and they can also be 
a candidate for the LKP |I2|. 

We are interested in investigating the inclusion of the graviton tower in the sce- 
nario with i?^ as the LKP. The gravitons do not evolve thermally, and have very long 
lifetimes, so if they are present, we expect them to decay to the LKP sometime after 
Big Bang Nucleosynthesis (BBN). The preservation of the light element abundances 
sets a bound on the amount of energy that can be released in a decay [131 CH CHI- 

^The radion can also be quite relevant ^Jj. We shall assume that the radion is inert during 
inflation and that it acquires a large mass, and therefore it does not have any impact on our 
analysis. 
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Secondly, since the gravitons decay ultimately into the LKP, we expect the relic 
density to increase. Therefore, the mass of the LKP consistent with the observed 
relic density is lowered. Additionally, if any of the gravitons decay after matter 
domination, we have to consider the effects of non-thermalized photons released in 
the decays on the spectrum of the diffuse photon flux (DPF) [16\ . 

This article is organized as follows. In Section 2, we review the standard calcula- 
tion for relic density. Section 3 reiterates the calculation for the density of gravitons 
presented in Ref. ^21; emphasizing the points relevant for our analysis. In Section 4, 
we analyze the lifetimes for the decay of the graviton tower and study its implication 
on the diffuse photon spectrum. Section 5 goes over the constraints on the energy 
released in the decay of the KK gravitons to the LKP. In Section 6, we combine the 
above results with the standard relic density calculation. Our goal is to analyze the 
effect of the gravitons on the predicted mass of the LKP consistent with the known 
dark matter density. We will show that, in fact, almost any value of lower 
than the one obtained in the absence of gravitons would be allowed, provided the 
reheating temperature, Tr, is large enough. For the calculation of the relic density 
of B^, we ignore complicating factors such as coannihilation and second KK level 
resonance effects jS[ El 13 El UH] ■ However, since these effects can be parameterized 
in the effective cross-section, we present the effect the graviton tower will have on 
the mass of the LKP without the gravitons, mwG, and note that as mwc increases, 
the contribution due to the graviton tower becomes large, and should be included in 
any precise calculation of mxK- We also find that all the gravitons, except for G^, 
decay right after BBN. At these early times, the electromagnetic BBN constraints 
are very weak. However it has been shown that at small lifetimes, the hadronic 
constraints are very stringent |17t I18j. Therefore, including constraints from both 
hadronic and electromagnetic decays, requiring consistency with BBN light element 
predictions, we find a stringent limit on the mass difference between mc-i and the 
LKP mass, mxK- Since the is long-lived, we also derive another constraint on 
the mass difference such that the observed spectrum of the DPF is not destroyed. 
Comparing the two constraints, we find that there exists a region of parameter space 
where both constraints are satisfied and which gives mass differences of the same 
order of magnitude as that obtained by radiative loop corrections to m^i jTH]. Ad- 
ditionally, we observe that if experimental evidence for UED is found, and the relic 
density induced by standard interactions of the KK modes is found to be lower than 
the one observed experimentally, the reheating temperature may be determined by 
assuming that the graviton KK modes provide the contribution necessary to achieve 
consistency between theory and experiment. Alternatively, including the possibil- 
ity of other unobserved exotic particles contributing to the relic density, an upper 
bound on the reheating temperature is obtained. 
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2 Relic Calculation 



Let us start by recalling the computation of the density of thermal relics. For any 
particle the evolution of the number density is governed by the Boltzmann 
equation: 



^+3Hn. = -{av){nl-n%) (1) 
H ^ i'-^f (2) 

- «(f )''»""'^ (3) 

where {av) is the thermally averaged cross-section multiplied by the relative velocity, 
and Eq. (0) is valid in the non-relativistic approximation, with m and g being 
the mass and number of degrees of freedom associated with the particle z. The 
temperature at which the particle decouples from the thermal bath is denoted by 
Tp {freeze-out temperature) and roughly corresponds to when T = n (av) is of the 
same order as H. 

Changing variables from n to Y = n/s (we will drop the subscript z from now 
on), where the entropy density is given by s = ^g^,T^, and using the fact that sR^ 
remains constant, we obtain. 



Ys = -{av)s\Y'-Y^\). (4) 



Introducing the variable x = the radiation dominated era: 



_ T^'g.T^ J_ dx_ 



s{Y'-Yl). (6) 



Then, Eq. (0} may be rewritten as 

dY_ _ {av) 
dx Hx 

We are interested in obtaining the rehc density, or the equivalent Foo, at late 
times. The solution of Eq. (jH)) enables the determination of the relic density and 
also the value of function of V. The ratio of the mass to the freeze-out 

temperature, xp = mi^i^/Tp is approximately given by 

g mKKMi{avy 



xp = In c(c + 2)v^^ -/J, ' (7) 
\ g^ X p / 
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and can be solved iteratively, with c ~ 1/2 and M4 = 1.7 x 10^^ TeV. As calculated 
in 0], for as the LKP, without including co- annihilation ^, the non-relativistic 
approximation is: 



a + 



6b 



X 



Anal {2Yf + 3Yb) 



9m^K 



KK 



I {2Yf + SFs) _ b' 



(9) 
(10) 



KK 



where Yf and Yb denote a summation over the fourth power of the hypercharges 
of fermions and bosons produced in the annihilation of respectively. Moreover, 
ignoring the splitting between the different KK states of a given tower, one obtains, 
approximately 

, [1.79 X 10^^ TeV / , 66'\1 
xp = In — a H . (11) 

niKK \ XfI 



1/2 

X p 

In terms of this, Yno is found to be: 



47r 



-1. 36 

^ M^rnKKXp (a H 

5 Xp 



3.05 X 10^6 TeV , , 36', 

(a' + — ) 

Xp rrixK Xp 



(12) 



where we have assumed masses of O (1) TeV, and used (7* ~ 92 for Tp ~ 50 GeV 
and So = 2889.2 cm~^. For more than one extra dimension this analysis is more 
complicated but a rough estimate may be obtained by assuming independent towers 
of KK modes for each extra dimension, d. The current relic energy density then is 
simply given by p = mxK Sq d Y^. We will however restrict ourselves to one extra 
dimension in most of our numerical work. 

As emphasized in the introduction, the above derivation for the density of B^ 
ignores the contribution of the KK modes of the gravitons, which may become 
relevant for sufficiently large values of the reheating temperature Tr. 



3 Density of KK Gravitons 

After reheating, except for gravitons, all other particles are initially in equilibrium 
and follow the pattern described in the previous section. After dilution of their den- 
sity due to inflation, the gravitons are produced mainly in collisions involving gauge 

^The effects of including these corrections to the cross-section in our analysis are discussed in 
detail in Section 6.3. 
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bosons, either in the initial or the final state, with cross-sections proportional to 
ai/Ml, where is the relevant gauge coupling. Since their density is so small, they 
are never in equilibrium and further we can ignore graviton annihilation processes 
to first order. Because the decay times for gravitons are of the order of r ~ 10'' s, 
the gravitons produced in the early universe are still present when the LKP freezes 
out. As the universe cools down further, the gravitons start to decay. Including 
higher order mass corrections pn); we will have decays of the form: 

n LKP + X (13) 

where X denotes SM particles. Here we are considering only KK number preserving 
decays since the KK number violating decays would only take place at the orbifold 
fixed points and so would only have a sizeable contribution in the case when the KK 
number preserving decays are suppressed due to phase space considerations. Hence, 
after decaying, the KK gravitons contribute to the LKP relic density and a proper 
computation of the dark matter density demands the addition of this effect to the 
relic density calculated in the previous section. 

In order to obtain a quantitative estimate of this effect, we follow the derivation 
for the number density of gravitons presented in Ref. [^. The number density of 
the KK gravitons at each level n is again determined by the Boltzmann equation: 

^ + 3HnG^ = CG^, (14) 

where Cq" is the collision operator and can be parameterized as follows: 

Cg^ = Ca[g4T)no]\ (15) 

where 

Here as is the strong coupling constant and C can be understood as the fraction of 
all possible collisions which will interact strongly to produce gravitons. The exact 
calculation of the total production cross section is quite complicated, but can be 
estimated by the method presented in Appendix |X| The authors of Ref. [12] used 
an analogy with the calculation of gravitino abundances [201 |2I] , and estimated C ~ 
O (0.01). However the sample calculation of the cross-section for the production of 
relativistic gravitons presented in Appendix pleads us instead to expect C ~ (9 (1). 
For completeness we will consider a range of values, 0.01 < C < 1, in our numerical 
work. 

For UED theories, g*(T) = g^^Dd(T), where g^^ is the effective number of 
degrees of freedom per KK level. Dd{T) can be approximated by counting all modes 
with masses below T: 



1 r 2^ nrf 



friKK 



(17) 
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where 

^^^ = 1X1:^ = 2'^' 3^' 2""'- ford = 1,2, 3..., (18) 

is the d dimensional volume of a unit sphere, and the factor 1/2'^ in Eq. (|17p accounts 
for the restriction to non-negative n. Assuming entropy conservation, S = sR^ oc 
g^{T)T^R^ oc T^+'^R^: 

Ids IdR dT 3 

-— = -3-— = -3if, — = -HT. 19 

sdt R dt ' dt 3 + d ^ ^ 

As in the previous section, making the substitution Y = n/s and with the help 
of Eq. ()19p the Boltzmann equation becomes: 

dYcr^ 3 + d 1 2 



dT 3 HTs 



Ca[g,iT)n,Y. (20) 



Yen changes until production stops at temperatures T ~ nrriKK and then 
remains constant until gravitons begin to decay. After BBN and before gravitons 
decay: 



45v^C'(3) rnxK 



27r8 



-"3- 



Ma 



3 + d Vd 
2 + rfV 2^ 



R 



rriKK/ 



i+- 



for n < 



R 



rriKK 



for n > 



R 



TrtKK 



(21) 



To find the Yq contribution to the LKP, assuming only KK number preserving 
decays. 



Yn 



A. 



Tr/tyikk 



nYc^ d'^n 



45v^C'(3) mKK 



27r« 



3 + d 



(x{d) C 
27r^/2 



TeV 



M4 • 



1 + c/)(4 + 3c/) V 23^^ 



mKK 



"^KK 



r(rf/2) 

where ~ 0.1 and gf^ ~ 200 



2,27r, ... for d = 1,2,... , 



(22) 
(23) 
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4 Decay Lifetimes and Diffuse Photon Flux 



The decay widths for the decay of the gravitons are calculated in Appendix El As 
shown there, the decay of G" for n > 1 is primarily into gauge bosons of KK number 
n/2 {{n ± l)/2 for n odd), assuming that all fermions are heavier. The lifetime for 
these is suppressed by powers of n, and is given by: 



327r 



rrir. 



v/2cos2 9w 1^1 V 



1.76 X 10^ s 



rTeVi 


5/2 


rTeVi 


1/2 


- rrin . 




[ An \ 





5.56 X 10^ s 



TeV 



(24) 



where A„ = mc" — 2m^^ << men and in the last line we have approximated 
An ~ nAi = n{mGi — itlbi) with Ai ~ lO^^rnxK; which is of the order of the mass 
corrections induced at one loop , and, as we will show below, is also of the order 
of the mass differences required to satisfy the phenomenological constraints coming 
from BBN and diffuse gamma ray constraints. 

The only long-lived graviton is G^, with a lifetime given by: 



371 Ml 
cos^ 6\Y A\ 

2.33 X 10^ s 
2.33 X 10^^ s 



TeV 



Ai 
TeV 



rrtKK 



(25) 



where we have again assumed in the last line that Ai ~ 10^ ttikk- The different 
dependence on the mass difference in Eq. and (^3)) is due to the fact that one 
of the decay products for is massless. 

The diffuse photon flux is only sensitive to photons released after matter dom- 
ination [THl 1221 since the late produced photons don't have time to thermalize. 
Therefore we just need to ensure that the energy released in the decay of the G^ 
does not destroy the spectrum. The diffuse photon flux is given by 



3c Ar^"i 



tn 



-7 J 



1/2 



3/2 



I0_ 



e(e^ - E, 



7/' 



(26) 



where is initial energy released in the decay, E^ is energy of photons observed 
now, to and Vq are the current time and volume of the universe, and A^^^i is the total 
number of G^s initially. Using Yqi calculated in the previous section: 



7 




Figure 1: The observed Diffuse Photon Fhix plotted in addition to the maximum flux from 
Eq. (|29|l for different values of C and tokk, assuming Tr such that f^si = 0.23. 



IK 



1.63 X MeV"^ 



Ai 



3/2 



.MeVj 

This has maximal value at: 



rriKK 



LTeV 



LMeV, 



1/2 



riKK/ 

f£2^y^'l.85xl0- 



e: 
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TGI \ 

3to. 



0.687 MeV 



GeV 



(27) 



(2^ 



Observe that, although is proportional to Ai, the inverse dependence of E"™"^ 
on Ai comes from the energy redshift induced by the expansion of the Universe 
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(tg-i oc '^). This, in turn, leads to a maximum value for the flux. 



3.06 X 10"^ sr"^ s"^ cm'^ MeV~^ 



3/2 



Ai 



GeV 



rriKK 



TeV 



• (29) 



If, for example, we assume a reheating temperature Tr ~ 50 mxK, we obtain 



dE.. 



1.1 X 10^2 sr^^ s"^ cm-2 MeV'^ 



Ai 



GeV 



^KK 

TeV 



(30) 



This is two orders of magnitude less than was obtained in Ref. ||^ for a similar mass 
range, and is marginally consistent with the observed diffuse photon flux ^^jj. The 
source of the difference is not computational, but based on the different assumptions 
made in both works. In Ref. j22|, the NLKP was assumed to be and had a density 
comparable to the dark matter relic density, while in our framework we are concerned 
with the decay of the G^s, which have a density much smaller than that associated 
with the dark matter relic density. 

In order to understand the constraints imposed by the photon flux on the 
decay, we can parameterize the experimentally observed flux in terms of E^ as 
follows. 



d^ 

dE^ 



4 X 10"^ sr"^ s"^ 



cm"^ MeV" 



MeV 



n 3 



E, 



7 J 



(31) 



This is plotted in Fig. Q For comparison we have also plotted the maximum flux 
according to Eq. for a range of values of C and mxK, assuming a Tr such that 
fi^i = 0.23. The mass difference allowed by the diffuse photon flux can then be 
quantifled by requiring that the maximum differential flux calculated from Eq. 
is less than the observed flux Eq. (jHH). It should also be noted that here we assumed 
that all the photons produced in the decay of the gravitons did not thermalize and 
are available to distort the spectrum now. This is certainly an overestimate. The 
constraints on Ai will be discussed in more detail in Section 6.6. 



5 Energy Released in Decays 

The electromagnetic showers produced by the decay of n LKP + X can create 

and destroy light elements. As discussed in [13 Cll EH] , this sets constraints on the 
lifetime r of the unstable as well as on the energy released per background 
photon in the decay. 

The dependence on r can be understood by looking at the characteristic energy 
scales in the initially produced photon spectrum. The primary photon created in 
the decay interacts with the background and creates an EM cascade. The fastest 
interactions are pair production and inverse compton scattering. These processes 
rapidly redistribute the energy and the non-thermal photons reach a quasi-static 
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equilibrium (QSE). The zeroth order QSE photon spectrum depends inversely on 
the temperature of the background plasma. If we make the uniform decay approx- 
imation, i.e. all particles decay at t = r, this corresponds to a cut-off energy of 
Ec = 103 MeV (r/10^ s)^/^ J^^. Therefore higher energy photo-erosion processes 
occur for longer lifetime values. 

We can understand the dependance of photodestruction and secondary pro- 
duction in a similar way. In the limit of small the decaying particle does not 
influence the light element abundances. Beyond this trivial case, we can again use 
the uniform decay approximation to gain some insights. As long as a reaction can 
take place {Eth, the threshold energy for a reaction, < Eq), a typical shower photon 
has energy: 

<^) = ^"'^^ (iow)"'(iiF7)"'- <^^) 

Therefore, the number of such photons per decay is ~ 5mKK/ (E), where Sttikk 
is the energy released in each decay. Thus, the non-thermal photon density is: 

= N^unlkp = '^,2 (33) 

These photons then further thermalize as well as cause photodestruction of par- 
ticles to yield other species of particles. Therefore the dependence of the fractional 
change in the abundance of a particle is proportional to and inversely proportional 
to {E) and therefore to r. 

Detailed numerical work gives a complicated but weak constraint on for getting 
the correct light element abundances for small lifetimes. However after about r > 
10'^ s, the dependance on steadies out to ^ < 10~^^ TeV [131 CH CS] • Comparing 
Eq. (J211) and (j2SI), we see that the only long-lived graviton is the with t{G^) ~ 



TeV 



, much larger than 10 s. Hence, electromagnetic energy release bounds 
only constrain the decay of the to the LKP. 

Similarly, hadrons produced in the decays can scatter off photons, electrons 
and, more importantly, background nuclei. Scattering off photons and electrons 
just causes them to lose energy, but the nuclei become more energetic and cause 
"hadronic showers" . Additionally if inelastic scattering occurs, the background nu- 
clei get dissociated and the light element abundances are changed. Due to the mul- 
titude of possible interactions, these are a lot more complicated to analyze than the 
purely electromagnetic case. Previously it had always been assumed that since the 
branching ratio into hadrons is small, the effect would also be negligible. However 
detailed numerical work f3CB] shows that in fact at early times, the constraints due 
to hadronic processes are much stronger then those due to electromagnetic processes. 
As can be seen from Eq. (j^^ the lifetime of the G"^s are suppressed by powers 

r n 3 

of n: t{G"') ~ 10^ s ^ . Therefore all the G"'s except for G^ decay right 
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after BBN. The bounds at these relatively early decay times depend strongly on the 
branching ratio of the decay of the gravitons into hadrons. If all the visible energy 
released is in the form of hadrons, the bounds are quite severe: C,h < 10~^^ TeV. 
If most of the decay is into leptons, then the bounds on C,h may be significantly 
weaker. A preferential decay into leptons could happen if, for any given tower, the 
KK leptons are very close in mass to the KK hypercharge gauge bosons, while the 
KK quarks are heavier (see Appendix El for calculation of the graviton decay widths 
into fermions). This would be the result one would obtain if including only the 
radiative corrections to the KK particle masses computed in Ref. ^Hl- Since this is 
very model dependent, we will again use the most constrained case of < 10~^^ 
TeV. 

Hence, both for early and late time decays we will use a conservative bound for 
the energy released per background photon, C,b < 10~^^ TeV, in our numerical work. 
Observe that while strong violations of this bound will certainly induce strong effects 
on the light element abundances, and therefore are ruled out, small violations of this 
bound are still possible due to the fact that we used very conservative limits on the 
energy released. 

One loop effects introduce corrections of two forms, constant (independent of 
n) and proportional to n (bulk and boundary correction terms) [TU] . Therefore, the 
mass of a particle at KK level n can be written as: 

71 

mn = -{l + S) + 6' (34) 
K 

where the corrections S and 6' for each particle type is the same at each level. The 
graviton KK modes, of course, do not receive corrections of this kind at the one loop 
level, and we shall assume, for simplicity, that their masses remain unperturbed: 
mQTi = n/R. 

To analyze the effects of the decays, we are going to treat the two cases (5, 6' = 0) 
separately. Note that this is only for ease of understanding. In fact, the true effect 
would just be the sum of the two terms. In the decay nB^ +X assuming that 

X is massless to first order, mgi ^ rriBi, and small momenta, the energy released 
in each decay is given by: 

^ - n^m|i ^ f nrriKKS for S' = , . 

^ 2mG" I for 5 = 0. ^ ' 

For this energy release to be consistent with the light element abundances ^3 
HH [T3] we need to have: 

e = / i?EM/Had^i?X^'^ <^B = 10"^^ TeV (36) 

Jo 

here nc^ = SoYcr^ as given in Eq. (j^U, and -BEM/Had is the electromagnetic/hadronic 
branching ratio. This fraction depends on whether the LKP is a lepton or a gauge 
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boson and the mass splittings at each KK level Throughout this article we 

choose -BEM/Had = 1 as the most stringent constraint possible. 

6 Determination of the Lightest KK Particle Mass 

In this section, we will discuss the determination of the lightest KK particle mass 
consistent with the observed dark matter relic density. Since the energy density 
of LKP is determined, in part, by the primordial graviton KK density, the lightest 
LKP mass mxK will depend on the reheating temperature and the parameter C 
governing graviton production at early times. 

6.1 Tr Reheating Temperature 

Including the decay of the gravitons in the density for B^, we have: 



Ybi = dY^ + Yc (37) 
ubi = SqYbi (38) 
Qb^ = (39) 

Pc 

where pc = 5.3 x 10^^ TeV/cm^, and Y^o is the abundance of the LKP without 
the inclusion of the gravitons, Eq. ()12|) . and Yq is the abundance of the gravitons, 
Eq. (|22p . We have not included either coannihilation or second KK resonance effects 
in the calculation of Y^o- 

Requiring that Qb^ — 0.23, we get : 



'r ^ rriKK 



(40) 



In principle, the above expression shows that the relic density constraint can be 
satisfied for any value of the lightest KK mode mass, rnxK; provided the value of 
Tr is adjusted according to Eq. (jlU)). In practice. Fig. |21 shows that in order to get 
a consistent dark matter density for low values of mxK, one would need large values 
of the ratio of the reheating temperature to ttt-kk, but as is discussed in the next 
section, these are disfavored due to strong coupling constraints. Similar results are 
presented in Fig. El for the case of D = 6. 

The mass of the LKP that would be obtained without the graviton tower con- 
tribution, mwG! is shown in the figure as the point where the bound on Tr becomes 
very strong, mxK ~ 0.92 TeV for D = 5 and mxK ~ 0.65 TeV for D = 6, since in 
this case the reheating temperature must become smaller than the lightest KK mass 
so that the relic density is not effected in any significant way by the gravitons. 
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C = 0.01 




"^KK (TeV) 



Fi gure 2: Values of the reheating temperature Tr obtained by demanding a proper dark matter 
density, ri^i = J^dm — 0.23, for different values of the graviton production parameter C = 1, 0.1, 
0.01, and assuming D = 5. Also shown are lines of constant ratios of the reheating temperature 
to the hghtest KK mass, Tr — 40 tokk, 70 niKK and 100 niKK- 



It is important to stress that for D = 5, values of the KK masses smaller than 
about 500 GeV, are also restricted in the minimal UED case by precision electroweak 
constraints [21 El 1211 ESj . Somewhat lower masses may be obtained for large values of 
the Higgs mass. However even then, within this minimal framework, values of mxK 
smaller than 350 GeV are strongly disfavored. From Fig. |21we see that for D = 5, 
C ~ 1 and values of the reheating temperature Tr < 40 mxK (consistent with 
the strong coupling bounds discussed in the next section), a consistent relic density 
may be obtained for any value of larger than 580 GeV. This range of values 
of the LKP mass are also consistent with precision electroweak constraints. Even 
for C ~ 0.1 and for the same range of values for the reheating temperature there 
could be modifications of up to ten percent of the value obtained without 
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Fi gure 3: Values of the reheating temperature Tn consistent with the requirement Clgi ~ f^DM = 
0.23, for C = 1, 0.1, 0.01, and assuming D = 6. Also shown in the Figure are curves of constant 
ratio of the reheating temperature to the lightest KK mode, Tn = 10 tokk, 20 ttikk and 30 tokk. 

the inclusion of gravitons. Larger (smaller) modifications are possible for larger 
(smaller) values of Tr. 

As shown in Fig. El for D = 6 one could obtain larger effects for smaller values 
of the ratio of Ts/rnKK- However, as will be shown in the next section, the bounds 
obtained from the requirement of perturbative consistency of the theory become 
much more stringent in this case. Therefore, large departures from the values 
obtained in the absence of gravitons seem to be disfavored for D = 6 within this 
minimal framework. 
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6.2 Constraints on the Reheating Temperature 

Large values of the reheating temperature compared to rriKK immediately raise the 
question of the ultraviolet cut-off for our effective 4D theory. In Ref . |2] , the authors 
assumed the limit on the KK masses to be the order of ~ 40 mxK, based on the 
quantum corrections to the strong gauge coupling. More stringent bounds have been 
evaluated in [211 12ZI I2H1 ■ Here we show that, if one computes the running of the zero 
mode hypercharge gauge coupling, and require it to remain weak in the ultraviolet 
regime, one obtains a bound similar to the one obtained in Ref. j2]. 

Let us then analyze the running of the zero mode gauge couplings for one extra 
dimension: 

In ^ - P^{Nkk HNkk] - HNkkI]) (41) 



9j ^(/^) 




where 






(81 ( 
1^10' 


6, 


= (-- 



57r^ ovr^ 



, , each KK level 

2 J ' 

SM (42) 

For D = 5, the most dangerous one for the KK scenario was the U(l) case and 
we found that gi develops a Landau pole at scales fi > 46 mxK [201 IHIIl IHI] • Our 
conclusion is then that above energy scales of about 40 mxK our effective theory 
breaks down. Similarly, for D = Q, a. Landau pole would develop at scales larger 
than about 10 mxK- However, these bounds may be avoided by assuming that 
this theory is just the low energy manifestation of a theory with more degrees of 
freedom and based on a higher, asymptotically free, gauge group. This would not 
only prevent ai from developing a Landau pole, but additionally it would increase 
the value of C. We have conservatively kept values up to Tr ~ 100 itikk for D = 5 
and Tr ~ 30 ttikk for D = 6. 

The impact of these bounds on the possible values of the lightest KK mass are 
displayed in Figs. El and El As we mentioned before, we see that if we restrict 
ourselves to the bounds implied by perturbative consistency, the mass of the LKP 
can be as low as ~ 0.6 TeV for both D = 5 and D = 6, but lower values, as low 
as those implied by consistency with precision electroweak measurements, would be 
only obtained by the possible relaxation of these bounds by new physics. On the 
other hand, if the most stringent bounds obtained in [211 12Z1 12H] were imposed, very 
small modifications to the mass of the LKP would be obtained. 
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Figure 4: Values of the ratio of the LKP mass consistent with JIdm ~ 0.23 to the one obtained 
in the absence of gravitons, mwG, for Tr — 20 rnxK, D = 5 and for C — 0.01, 0.1 and 1. 

6.3 Additional Contributions to the Annihilation Cross Sec- 
tion 

Until now, we have only considered the effects induced by the annihilation of -B^s 
on the total relic density. However, several effects can change the density of the 
B^s. For example, if we assume that the mass difference between the LKP and the 
NLKP is less than about 10%, the NLKP is still thermally accessible to the LKP, and 
the effective cross-section, including coannihilations must be used in Eq. (Q). Also, 
possible resonant effects induced by the second KK-level modes must be included 
13 IHl EDI- Note that even though including these effects changes the mass of the 
LKP consistent with the dark matter density without the inclusion of the gravitons, 
mwG, and will certainly change any numerical results we obtain, the qualitative 
picture we have presented here will remain unchanged. 
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Figure 5: Same as Fig. |3| but for Tr = 40 ttikk- 

The corrections to the cross-section can be quantified by considering a' and b' in 
Eq. (fT^ to encode all the information about the processes that affect Yoo- Due to 
the weak, logarithmic dependence of on mxK and the interaction cross section, 
we can approximately parameterize Yoo as being proportional to rnxK, 



y 



/mKK\ 
V TeV y ■ 



(43) 



Then y, such that mwG (including coannihilation, second KK resonances etc., but 
without the graviton tower) is consistent with the experimentally observed dark 
matter density, floM, is given by: 



y 



DM 



(44) 



We can now replace Fqo using Eq. (jlHI) and (jl^ in Eq. (j?fj) . Using this, we get 
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Fi gure 6: Same as Fig. 0] but for Tr — 100 tokk • 



the change in mwc necessary to reproduce the observed dark matter density for a 
given reheating temperature Tr: 



^DMPc 



1/2 



niKK 



4+3d 
2 



(45) 



This is plotted in Fig. |1[ Elandinifor Tr = 20 rnxK, 40 mxK and 100 mxK and 
a range of mwc predicted in Refs. El 13 El HD]- We see that even for Tr = 
20 rriKK with C = 1, the mass consistent with the observed rehc density undergoes 
a relevant modification due to the inclusion of the gravitons. This effect becomes 
increasingly important with increasing m-wc- Therefore, depending on the reheating 
temperature, any precise calculation of niust include the contribution from the 
graviton tower to be accurate. 



18 



It is interesting to observe that for Tr ~ 100 ttikk and C ~ 1, the necessary 
values of mxK are below the bounds imposed by precision electroweak constraints. 
Conversely, this shows that for this value of C and for values of the KK masses 
consistent with precision electroweak constraints, such large values of the reheating 
temperature will lead to an excess of dark matter density and are therefore disfa- 
vored. On the contrary, for values of Tr < 40mKK, the obtained values of m-KK are 
in good agreement with the scale set by precision electroweak data. This shows an 
interesting correlation between the dark matter density bounds and those coming 
from requiring the perturbative consistency of the theory up to scales of the order 
ofT^. 



6.4 Constraints on the G —B Mass Difference for 5' = 

In the discussion above, we have ignored the bounds on the energy released in the 
graviton decays. As we have stated in section 4, the amount of energy released will 
depend on the type of mass correction the graviton and receive (see Eq. (jH^j) 
and Eq. (jHBj) ). For the case of 5' = 0, the energy released per background photon is 
given by 

e = ^^KK ^0 Ya<iB = 10-^^ TeV. (46) 

Assuming, as stated before, that there are no corrections to the KK graviton masses 
mG<7i, then = nniQi . In addition we can identify with m^i . The corrections 
to m^i = mQi{l — 6) are then equivalent to corrections to men = nmBi(l + 6) for 
small 6. 

If we now assume that the reheating temperature is such that one reproduces the 
correct dark matter relic density, ^2^1= 0.23, as obtained in the previous section, 
then we get a constraint on the mass difference between and B^: 

This is plotted in Fig. [3 for D = 5. This constraint is independent of C. The mass 
difference between B^ and is very tightly constrained to be less than a GeV 
for most of the mass range considered. As approaches the value of the mass 
consistent with a proper dark matter relic density in the absence of gravitons, the 
bounds become very weak. This is due to the fact that in this mentioned 
before, Tr may be of the order of (or smaller than) mxK, and hence very few 
gravitons would be produced, implying that the mass difference can be large without 
inducing any dramatic effects. 

For most of the parameter space, we have derived here a very stringent constraint 
on the mass difference. Therefore, we need to consider whether such small mass 
differences would require too much fine tuning to be considered natural. As is 
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Figure 7: Values of the allowed mass difference between and the LKP, Sttlkk, for S' — 0, 
assuming f^oM = 0.23 and ^ < 10^^^ TeV for D = 5. We have ignored possible coannihilation and 
second KK level effects. 



discussed in Section 6.7, we find that in fact the mass differences obtained here are 
the order of magnitude of the one loop corrections induced to the mass. 

6.5 ^ = 

Here we depart from the previous assumption and assume that all graviton KK 
modes receive a constant, positive mass correction with respect to nrriKK- The 
energy released for this kind of mass correction tends to be much smaller than in 
the case of 5 7^ analyzed above (see Eq. (jlUj) ). and is given by: 



B 



EM/Had 



S'soY^ <^B = 10"^^ TeV 



(4J 
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Figure 8: Values of 5' if 17dm = 0.23 and ^ < TeV for Z? = 5, C = 1, 0.1 and 0.01. We 

have ignored possible coannihilation and second KK level effects. 



= J Ycnd'^n 



VdA^, 3 + d 



27r8 



P{d) C 



23*^ d{2 + 3d) 



R 



rriKK 



1+' 



TeV 



R 



mKK 



'-^ 2 



(49) 



For this case, since the dependance of Yq on Tr is different from that of Q, we get 
a more comphcated expression for 6', depending exphcitly on the number of extra 
dimensions and C: 



6' < 



BEM/HaxiPid)C SomKK 



a{d)soCm'^^ 



^ 2+3d 
4+3d 



(50) 
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This is plotted in Fig. |H1 As anticipated, the bounds on 5' are much weaker than 
on (5mKK- However, even though the constraint is much weaker now than before, 
only small mass differences, smaller than about 10 GeV, are allowed in most of the 
parameter space. 



6.6 Diffuse Photon Flux Constraints on the 
Difference 



Mass 



Since the decays late to the we have to consider its effect on the diffuse photon 
flux. As long as the decays after matter domination, it can have an effect on 
the diffuse photon flux. Since the observed flux is smooth to a high degree, we can 
derive a constraint on the mass difference Ai between and by requiring that 
Eq. be less than Eq. (j^ . Using this inequality, we are led to: 



Ai > 2.48 X 10"^ 



T, 



R 



rriKK 



3/2 



1 



rnxK 



TeV 



1/2 



GeV, 



(51) 



where Tr is given by Eq. This is plotted in Fig. El We see again that as wkk 

approaches the value consistent with the dark matter density without gravitons, the 
constraint on Ai becomes very weak. At these masses, even if the were decaying 
right now, their density is so small that the decays are basically invisible and don't 
affect the photon spectrum. 



6.7 Comparison of Constraints 

Comparing the constraints on the mass difference between B^ and G^, plotted in 
Fig. EH we see that there is a region of parameter space where both conditions are 
satisfied. The allowed region increases for larger values of G, and for G = 1, extends 
up to values of order 600 GeV, which is about the value of the bound placed by 
precision electroweak tests for a light Higgs boson. The range of allowed values of 
Ai shrinks for decreasing values of rnxx/'^WG- For instance, in the case represented 
in Fig.^Jfor G = 1, values of — 0.8 TeV may be obtained for a range of values 
of Ai ~ O (1) GeV, while values as low as mxK ~ 0.6 TeV are only allowed for 
Ai ~ 0.6 GeV. Interestingly enough the allowed values of ttikk agree with the ones 
consistent with the observed dark matter relic density for Tr < 40 rriKK depicted in 
Fig. El 

The effects on the allowed values of Ai as mwc changes are illustrated in Figs.lTD 
and^l As shown in these figures, the qualitative behavior for the different mwc is 
the same as the one found in Fig. ^1 
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Figure 9: Constraints on the value of the mass difference between and the LKP, Ai, assuming 
f^DM = 0.23 ioi D — 5, C = 1, 0.1 and 0.01. We have ignored possible coannihilation and second 
KK level effects. 



6.8 One Loop Corrections to the KK Masses 

We want to compare the bounds obtained on Ai due to BBN and the diffuse photon 
ffux found in the previous section, with the one-loop corrections induced in the 
mass, assuming that m^i = 1/R. Radiative corrections induced by boundary terms 
at the orbifold fixed points affect the spectrum of the standard KK particles. Since 
radiative corrections are divergent they may be regularized by counterterms localized 
at these boundaries. The possibility of including localized counterterms also implies 
that one could have started with different values of the localized boundary terms 
than the ones assumed in the minimal scenario, and this would lead to modifications 
of the spectrum different from the one discussed below [22] ■ We include, however, 
a discussion of this minimal scenario in order to show that, interestingly enough, 
the magnitude of the one-loop induced corrections is indeed of the order of the one 
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Figure 10: Constraints on the value of the mass difference between and the LKP, Ai, 
assuming fioM = 0.23 for D = 5, C — 1, 0.1 and 0.01. We have ignored possible coannihilation 
and second KK level effects. 



necessary to satisfy the BBN and diffuse gamma ray bounds. 

In order to compute the LKP mass, apart from the one-loop corrections, one 
must also take into account the gfv'^/A correction to the gauge boson mass matrix 
due to the vev of the Higgs field [12]. In the 5", basis: 



V 



,-^22 



\ 



9i92V 



4 



9i92V 



n 



+ 5{m 



2 ^ 



I 2 2 

+ -^92V 



(52) 



Since the mixing is very small compared to the mass difference between the 
and the induced at the one-loop level, the neutral gauge bosons approximately 
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Figure 11: Same as Fig.[Tni but for a value of mwG = 1.5 TeV. 



become pure and W^. Therefore, the mass correction to is very well approx- 
imated by: 



5(m 



39aiC(3) lai 



+ 



9W 



D=5 



(53) 



2 47r3 ' 6 47r /i" 

This correction to becomes positive for masses below about 800 GeV, and 
the absolute value is plotted in Fig. UHl Therefore, as noted in ilQ^, below this mass, 
assuming only one-loop corrections, the graviton is the LKP. This case has been 
extensively studied and the constraints are derived in Ref. |^. The S^s, with a 
density comparable to O ~ 0.23, would be decaying late to and based on an 
analysis similar to the one performed in Section 4, one concludes that, apart from 
the finely tuned case in which the mass difference approximately vanishes, small 
mass differences between the B^ and the would lead to a large impact on the 
diffuse photon flux. Therefore, unless some evidence is seen in the diffuse photon 
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Figure 12: Same as Fig. [113 but for a value of mwG = 2 TeV. 

flux for new pliysics for values of i?^ ~ (9 (1) MeV (see Fig.Q), this minimal scenario 
would be ruled out for mxK < 800 GeV ^. 

We see that in the allowed region demarcated in Figs. [TUIIIII and [T21 the values of 
Ai are, as emphasized before, of the same order of magnitude as the range predicted 
by the one loop corrections shown in Fig. ^1 The precise quantitative constraints 
on the possible values of mKK/''^WG depend on mwc- While for mwG — 1, 2 TeV, 
only corrections of the order of ten percent would be allowed, for mwG — 1-5 TeV, 
the one-loop corrections are consistent with those necessary to satisfy both the BBN 
and diffuse photon flux constraints ^. 

In evaluating these constraints we have to stress again that we have used a very 

^Please see Ref. [SJ for a possible solution to this problem, by the introduction of Dirac neu- 
trinos. 

■^In Ref. jlO| it was argued that, within this minimal framework, values of mwG > 1-4 TeV 
would be disfavored since a charged Higgs would become the LKP. 
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Figure 13: One loop corrections to the mass as a function of ttikk- 

conservative estimate for the hadronic branching ratio of Bh — 1. The actual energy 
released into hadrons is approximately proportional to Bh and would be somewhat 
lower since the KK quarks are much heavier then the KK right handed leptons. 
Additionally there are large errors in the detection of the diffuse photon flux. A 
proper computation of the hadronic branching ratio of the decay of gravitons, as 
well as a more accurate spectrum for the diffuse photon flux, would be necessary 
in order to determine the compatibility of this minimal scenario with the energy 
release constraints. 

6.9 Determination of the Reheating Temperature 

Let us finish this section by commenting on the possible experimental probes of 
this scenario. Due to their small couplings to matter, provided is not the LKP, 
the KK gravitons will not be produced at laboratory experiments. Therefore, their 
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Figure 14: The reheating temperature, Tr, required to make up the deficit dark matter density 
by KK gravitons as dictated by Eq. (|54|l . 



presence may only be probed by indirect effects, lilce their impact on cosmology as we 
have discussed in this article. In order to be able to evaluate these effects, one would 
have to find conclusive evidence of UED at, for instance, collider experiments and, 
in addition, measure the properties of the relevant first and second KK modes that 
contribute to the annihilation cross section. Then, under the assumption that 
no other exotic particle contributes to the dark matter density one could estimate 
the relic density associated with the LKP for the specific value of ttikk measured. 
If the obtained relic density disagrees with the experimentally measured one, then 
one would get information about the possible KK graviton contribution. 

Knowledge of the graviton contribution and the KK mass mKK will, in turn, 
allow us to determine the reheating temperature. By analogy with Eq. ()40j). the 
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reheating temperature for one extra dimension would be given by, 



( 



) 



2/7 



Tr = niKK 



a(l)Csom^K 



(54) 



where Af2 is the difference between the measured value of Q and the theoretically 
estimated one from the measured value of ttikk and the corresponding annihilation 
cross section. This is plotted for C = 1 in Fig. for several different values of AQ 
for reference. Even including the possibility of exotic particles contributing to the 
dark matter density, Eq. (j54j) gives an upper limit on the reheating temperature. 

7 Conclusions 

In this article, we analyzed the effects of including the KK graviton tower on the 
determination of the relic density associated with the lightest KK particle in the 
scenario of Universal Extra Dimensions. Gravitons may be copiously produced in the 
early universe and their subsequent decay into the LKP may lead to a modification 
of the observed LKP relic density. Graviton production is governed by a parameter 
C ~ and by the reheating temperature, which is bounded from above by the 
requirement of maintaining the perturbative consistency of the theory, Tr < 40mKK 
(Tr < lOrnKx) for D = 5 {D = 6). In the case that Universal Extra Dimensions are 
observed in laboratory experiments, we show that an upper limit on the reheating 
temperature can be deduced from the requirement that only the LKP contributes 
to the observed dark matter density. 

Throughout this work we have assumed KK parity conservation, leading to the 
stability of the LKP particle. We found that including the graviton effectively 
lowers the lightest KK particle mass consistent with the observed dark matter relic 
density from values of about 0.92 TeV to values as low as ttikk ~ 0.58 TeV for 
C = 1, D = 5 and Tr < 40 ttikk- Additionally, including effects which change the 
LKP density (coannihilation, second KK- level mode resonant contributions, etc.), 
we show that the graviton tower has a large impact on the predicted mass of the LKP 
with increasing mass. It should be stressed here that these results are independent 
of which KK particle is the LKP (or the NLKP). 

Additionally, there are bounds on the mass difference of the LKP with the gravi- 
ton KK modes induced by the requirement that the energy released in the graviton 
decay does not lead to a disturbance of the light element abundance or the diffuse 
photon flux. Under the assumption that the graviton spectrum is ~ n/R, we 
have obtained a bound on the mass difference of and the LKP (B^) mass, that 
is consistent with the minimal one-loop corrections obtained in Ref. ^H] for a large 
range of values of mxK- 
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(P^) 





Figure A. 1: Feynman Diagrams contributing to q''+q^ — > g^+G". The corresponding amplitudes 
are given m l|X5|) and lfO)l . 

APPENDIX 



A Graviton Production Cross-Section 

We will calculate explicitly the graviton production cross-section for q\ki)+q^{k2) 
g^{pi) + G^{p2). The feynman rules we will need are jl2[ I33j: 



X 



Sn,\l-m\ — 7^Sn,l+r. 



(A.l) 



g^ih) + G 



2Ma 



- [r]^'^k\k'i + r^'^^fcf - t^k^ikl^ + (z/ ^ /i) 



5. 



n,\l—m\ 



(A.2) 
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q(kj 



G" (Pj, 




g-'^" (k+k^) 



Figure A. 2: Feynman Diagrams contributing to q^+(f"- g^+G"^ . The corresponding amplitudes 
are given in l|A.7|l and l|A.8|l . 



q + q g: 



(A.3) 



igt"- 
AM. 



(A.4) 



Four diagrams contribute to this process (Fig. lA.lf IA.2jl . The amphtudes are 
i;iven by: 



igt" 
2Ma 



(A.5) 



4M 



2{ki.p2 - mirrir, 



(A.6) 
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(A.7) 



: 



(Ai 



In the center-of-mass frame, the cross-section is given by 



CM 



1 


P2 






yi - V2 


{2'k)HEcm 



\M{h,k2-^Pi,P2)\. (A.9) 




Here, Ei, E2 and z/i, 1/2 label the initial particle energies and velocities, and p2 is the 
three-momentum of the graviton. Since we are interested in the relativistic limit, 
we can use the more general form for all four particles with m = 0: 

Averaging over initial spins and colors and summing over final polarizations, \}A\^ 
was calculated using the Tensorial 3.0 package for Mathematica jSSl- The cross- 
section is then given by: 

«3 s 1184-M057r / a -, -, x 

< au >= %.d(G)C(r)7r . (A.ll) 

47rM| ^ ' ^ ' 1296 ^ ' 

The fraction of strongly interacting degrees of freedom per KK level is computed 
by counting the different allowed interactions. Four strong interactions can produce 
gravitons: 

q + q^ 9 + G 

g + q^ q + G 
g + q-^q + G 
9 + 9 ^ 9 + G 

There are 2 x 3 x 24 degrees of freedom for the quarks and 3 x 8 for the gluons. The 
fraction can then be calculated: 



p,2 = 24 



7\V 6 \\ ^/7\ f 2A \ f 6 



8/ V197.5/ V8/ V197.5/ V197.5 



/ 24 \ ^ 
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Assuming that all the other graviton production processes have a cross-section 
similar to ()A.11|) . C used to parameterize the graviton production cross-section is 
given by: 



^ ^ 1184+ IOStt ^ 
C = ps2d{G)C{r)'K = 2.79, 



1296 



(A. 13) 



B Decay Lifetimes 

The graviton polarization sum is given by [T^ : 



B 



mt 



mt 



4 
3 



kpkfj 



mt 



(B.i: 



The decay width of G^{p) f\ki) + /™(A;2) summed over final and averaged 
over initial polarizations is: 



I I fra\ 



|ki| 1 



Mml 5 ^, 

1 "^n 

1607r Ml 



-Y,u\k2)Y^^-v'' {k{)v'' {k{)Y*^f'u\k2)B^ 



fj,uaf3 



mt 



^ _ (mi - m^Y ^ _ {nil + ^m) 



mt 



5/2 



^ _^ 2 {mi - mn,)^ 



(B.2) 



This is suppressed for / ^ m. Therefore assuming that the main decay is into 
I = m = n/2: 



m=n/2 1 £'m=n/2\ 



1 "^n 

1607r M| 



4m? 



mt 



2 i5/2 



(B.3) 



The decay width of G"'{p) B\ki) + B'^{k2) summed over final and averaged 
over initial polarizations is given by: 



kll l^^^a/3^*paA« /^^^ 



^vrm^ 5 



^^2k 

m 



2 ; -BlJ'J^pcr 
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cos^ 9w rnl^ 



X 



SOtt 
13 



Mi 



1 - 



{nil 



+ 



7 {nil 



mt 



rrir. 



5/2 



1 - 



{nil 



nir 



1/2 



+ 



mt 



1 {mi + mmf 



ml 



(B.4) 



We see again that the decay into KK modes other than n/2 is suppressed. Therefore 
for all n > 1, assuming now that the main decay is in fact into 5"/^ + i?"'/^ with an 
extra factor of 1/2 for identical particles in the final state: 



with A„ = m„ — 2mm ~ "t.g" ~ "^B" << ^n- 



cos^ 9w 


ml 


\ 4m2; 


1/2 - 


SOtt 


Ml 


ml _ 





13 14m^ m"^ 



12 



3 ml 



v^cos^ 9w ml fK^ 
327r Ml^mn 

For any n decay into a KK gauge boson and a photon: 

cos^ 9w r„ ^m 



307r Ml 



6 + 3- 



m 



2 



+ 



m; 



mi 



(B.5) 



(B.6) 



(B.7) 
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